Magnetic Ordering of Nuclear Spins in an Interacting 2D Electron Gas 
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We investigate the magnetic behavior of nuclear spins embedded in a 2D interacting electron gas 
using a Kondo lattice model description. We derive an effective magnetic Hamiltonian for the nuclear 
spins which is of the RKKY type and where the interactions between the nuclear spins are strongly 
modified by the electron-electron interactions. We show that the nuclear magnetic ordering at finite 
temperature relies on the (anomalous) behavior of the 2D static electron spin susceptibility, and 
thus provides a connection between low-dimensional magnetism and non-analyticities in interacting 
2D electron systems. Using various perturbative and non-perturbative approximation schemes in 
order to establish the general shape of the electron spin susceptibility as function of its wave vector, 
we show that the nuclear spins locally order ferromagnetically, and that this ordering can become 
global in certain regimes of interest. We demonstrate that the associated Curie temperature for the 
nuclear system increases with the electron-electron interactions up to the millikelvin range. 

PACS numbers: 71.10.Ay,71.10.Ca,71.70.Gm 



I. INTRODUCTION 

In the last decade, the field of spintronics has seen re- 
markable progress pii^i^ Among them, the possibility of 
confining electron spins in quantum dots opens the door 
to quantum spintronics. This is based on the possibil- 
ity of controlling and manipulating single electron spins 
in order to build devices able to achieve operations for 
quantum information processing. The most promising 
and challenging idea is the use of spins of confined elec- 
trons in quantum dots to realize quantum bitsj^ Within 
the last years, all the necessary requirements for spin- 
based quantum computation have been realized experi- 
mentally, going from the coherent exchange of two elec- 
tron spins in a double dot^ to the coherent control of 
a single electron spin, including the observation of Rabi 
oscillations.— These achievements have become possible 
because electron spins in semiconductor quantum dots 
are relatively weakly coupled to their environment and 
therefore long lived quantities, quite robust against de- 
cay. Indeed, longitudinal relaxation times in these sys- 
tems have been measured to be of the order of Iseci^'^'^ 
A lower bound on the spin decoherence time for an en- 
semble of electron spins in GaAs quantum dots has been 
measured to be typical larger than 100 ns,— while a co- 
herence time in a single quantum dot exceeding 1 fis has 
been recently achieved using spin-echo techniquesi^ It is 
by now well established that one of the major sources 
of decoherence for a single electron spin confined in a 
quantum dot is the contact hyperfine interaction with 
the surrounding lattice nuclear spins.— 

One possibility to lift this source of decoherence is the 
development of quantum control techniques which effec- 
tively lessen or even suppress the nuclear spin coupling 
to the electron spin.— '^^'^^ Another possibility is to nar- 
row the nuclear spin distribution,— li^ii^ or dynamically 
polarize the nuclear spins i^^i^"^:^^:^^:^^ However, in order 



to extend the spin decay time by one order of magnitude 
through polarization of the nuclear spins, a polarization 
of above 99% is required,— quite far from the best result 
so far reached in quantum dots, which is around 60%.— 
A common point to the aforementioned approaches is 
their aim at mitigating nuclear spin fluctuations by ex- 
ternal actions. Recently, the possibility was raised of 
an intrinsic polarization of nuclear spins at finite but low 
temperature in the two dimensional electron gas (2DEG) 
confined by the GaAs heterostructurei^ 

The nuclear spins within the 2DEG interact mainly 
via the Rudermann-Kittel-Kasuya-Yosida (RKKY) 
interaction;^ which is mediated by the conduction 
electrons (the direct dipolar interactions between the 
nuclear spins are much weaker, see below). An intrinsic 
nuclear spin polarization relies on the existence of a 
temperature dependent magnetic phase transition, at 
which a ferromagnetic ordering sets in, thus defining a 
nuclear spin Curie temperature. 

The first estimate of such a Curie temperature was ob- 
tained for three-dimensional (3D) metallic samples, using 
a Weiss mean field treatment by Frohlich and Nabarro 
more than sixty years ago*^ They determined the nu- 
clear spin Curie temperature to be in the microkelvin 
range or less for 3D metals. A Weiss mean field treat- 
ment also gives a nuclear spin Curie temperature Tc in 
the microkelvin range for a typical 2DEG made from 
GaAs heterostructures,^° yet a more detailed analysis 
is desirable for at least two reasons. First, a Weiss 
mean field analysis does not take into account properly 
the dimensionality of the system, and second ignores 
electron-electron (e-e) interactions. In two dimensions 
(2D), the Mermin- Wagner theoremS2, states that there is 
no phase transition at finite temperature for spin sys- 
tems with Heisenberg (isotropic) interactions, provided 
that the interactions are short-ranged enough. However, 
RKKY interactions are long-ranged and strictly speak- 
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ing, the Mermin- Wagner theorem does not apply, al- 
though a conjecture extending the Mermin- Wagner the- 
orem for RKKY interactions due to non-interacting elec- 
tron systems has been recently formulated (and proved 
in some particular cases) 

In Ref. [IS, we started from a Kondo lattice descrip- 
tion for the system composed of nuclear spins and elec- 
trons, then derived a rather general effective Hamiltonian 
for nuclear spins after integrating out electron degrees 
of freedom, and finally performed a spin wave analysis 
around a ferromagnetic ground state (which we assumed 
to be the lowest energy state). We indeed showed that 
Tc = for non-interacting electrons in agreement with 
the latter conjecture. However, taking into account e-e 
interactions changed drastically this conclusion. It turns 
out that e-e interactions modify the long range nature 
of the 2D RKKY interactions (which are directly re- 
lated here to the static electron spin susceptibility) and 
thereby allow some ordering of the nuclear spins at finite 
temperaturei^ Furthermore, we showed that the temper- 
ature scale at which this ordering takes place is enhanced 
by e-e interactions .~ 

The study of thermodynamic quantities in interacting 
electron liquids (especially in 2D) has attracted some 
theoretical2^>2^>21i2ai2ai^i2ii^iMiM and experimental^ 
interest recently with the goal to find deviations from the 
standard Landau-Fermi liquid behavior. It is therefore 
quite remarkable that the macroscopic magnetic proper- 
ties of nuclear spins in a 2DEG, and thus their finite tem- 
perature ordering, are directly related to the corrections 
to the the static electron spin susceptibility induced by e- 
e interactions. They may therefore be associated with an 
indirect signature of Fermi Hquid non-analyticities. Nev- 
ertheless, it turns out that the temperature dependence 
of the electron spin susceptibility XsiT) is rather intri- 
cate. On the one hand, from perturbative calculations in 
second order in the short-ranged interaction strength one 
obtains that |xs(T)| increases with temperature . 
The same behavior is reproduced by effective supersym- 
metric theories.— On the other hand, non-perturbative 
calculations, taking into account renormaHzation effects, 
found that |xs(T)| has a non-monotonic behavior and 
first decreases with temperaturei^i^ This latter behav- 
ior is in agreement with recent experiments on 2DEGs.— 

In view of these recent controversial results, we want to 
reconsider the question of a finite temperature ordering 
of nuclear spins by taking into account renormaHzation 
effects of the static spin susceptibility Xs{q), where q is 
the wave vector, and therefore going beyond Ref. H^. It 
turns out that, a priori, different nuclear spin orderings 
can occur, depending on temperature and other sample 
parameters such as the interaction strength, measured 
by the dimensionless parameter Ts (essentially, the ratio 
between Coulomb and kinetic energy of the electrons). 
We consider at least two possible ordered phases in the 
nuclear system: a ferromagnetic orderingSi but also a 
helical spin ordering where the nuclear spins align ferro- 
magnetically at the scale of the nuclear lattice constant 



but point in opposite directions at the scale of the Fermi 
wave length (roughly two orders of magnitude larger than 
the nuclear lattice spacing at small r-g < 1). Depending 
on the general non-perturbative shape of Xsiq) (which 
may have a complex dependence on T and r^), we discuss 
the possible ordered phases and their associated magnetic 
properties. 

The outline of the paper is as follows: In Sec. m we for- 
mulate a Kondo lattice description of our problem where 
the nuclear spins are playing a role analogous to magnetic 
impurities embedded in an electron liquid. We then de- 
rive a general effective magnetic Hamiltonian for nuclear 
spins where the interaction is controlled by the electron 
spin susceptibility in 2D. In Sec. IIIIl we calculate the 
electron spin susceptibility in an interacting 2DEG using 
various approximation schemes for both, short-ranged 
and long-ranged interactions. Particular attention is paid 
to renormaHzation effects in the Cooper channel which 
turn out to be important. Sec. [IV] is devoted to the 
magnetic properties of the nuclear spins depending on 
the general wave-vector dependence of the electron spin 
susceptibility. We discuss two different phases: A ferro- 
magnetic phase and a heHcal phase with a period of the 
order of the electron Fermi wavelength. Finally, Sec. |V] 
contains a summary of our main results and also perspec- 
tives. Appendix[A]contains some details of the derivation 
of the effective nuclear spin Hamiltonian and of the re- 
duction to a strictly 2D problem. 



II. MODEL HAMILTONIAN 

A. Kondo lattice description 

In order to study an interacting electron gas coupled to 
nuclear spins within the 2DEG, we adopt a tight-binding 
representation in which each lattice site contains a single 
nuclear spin and electrons can hop between neighboring 
sites. A general Hamiltonian describing such a system 
reads 

1 ^' 

= Hq + Hn + Hdd, (1) 

where Hq denotes the conduction electron Hamiltonian, 
Hn the electron-nuclear spin hyperfine interaction and 
Hdd the general dipolar interaction between the nuclear 
spins. Hq can be rather general and includes electron- 
electron (e-e) interactions. In Eq. |[T|, c]^ creates an 
electron at the lattice site with spin a =t,i, and 
T represents the PauH matrices. We have also intro- 
duced Ij = the nuclear spin located at the 
lattice site r^, and Aj the hyperfine coupling constant 
between the electron and the nuclear spin at site Vj. 
Summation over the spin components a, (3 = x,y,z is 
implied. The electron spin operator at site is defined 
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by Sj = hcj^Tca'Cja' (for convenience we normalize the 
spin operator here to 1). Ni denotes the total number of 
lattice sites. >From here on, we assume that Aj — A > 0, 
which means we assume the hyperfine interaction is anti- 
ferromagnetic and the same for all atoms that constitute 
the heterostructures (typically Ga and As and their iso- 
topes). 



The nuclear spins are also coupled via the dipolar in- 
teraction to other nuclear spins, which are not embedded 
in the 2DEG. Taking into account this interaction as well 
makes the problem of the magnetism of nuclear spins in 
GaAs heterostructures an a priori 3D tremendously com- 
plicated one. Nevertheless, it turns out that the dipolar 
interaction energy scale Em is the smallest one. It has 
been estimated to be E^d ~ 100 nKM In particular, 
ksT ^ Edd, where T is the temperature of a typical ex- 
periment. In the rest of the paper, we neglect all direct 
dipolar interactions between the nuclear spins, which are 
in general smaller than the indirect interaction, as we will 
see. Therefore, we assume that vff « in Eq. ^. This 
assumption is important since it allows us to focus only 
on those nuclear spins which lie within the support of the 
electron envelope wave function (in growth direction). 

The general Hamiltonian in Eq. |(T]) is the well-known 
Kondo lattice Hamiltonian (KLH), though Hq contains 
also e-e interactions. The KLH is one of the most stud- 
ied models in condensed matter theory due to its large 
variety of applications. The KLH has been used to de- 
scribe the properties of transition metal oxides,— heavy 
fermions compounds,—"^ more recently also magnetic 
semiconductors (or semi-metals) in the series of rare 
earth substances;^, and diluted magnetic semiconduc- 
tors such as Gai_xMnxAs,— to list only a few. The 
nuclear spins play a role analogous to magnetic impuri- 
ties in the Kondo lattice problem. The regime in which 
we are interested corresponds to the weak Kondo cou- 
pling regime in the sense that A <^ Ep, where Ep is 
the Fermi energy. Furthermore, the nuclear spin den- 
sity ris is far larger than the electron density n^. It is 
worth noticing that the single nuclear spin Kondo tem- 
perature Tk ~ Dexp{—EF/A) (with D being the elec- 
tron bandwidth) is extremely small compared to all other 
energy scales. We are therefore far away from the so- 
called controversial exhaustion regime^ where the indi- 
vidual screening of the impurity competes with indirect 
magnetic exchange between the nuclear spins. 



In this low electron density regime, the ground state 
of the magnetic system (here the nuclear spins) has been 
shown to be ordered ferromagnetically in 3D using var- 
ious treatments that go beyond mean field theory and 
which notably include spin wave modes (but neglect e-e 
interactions)!^ 



B. Derivation of an effective magnetic Hamiltonian 

We first go to Fourier space and rewrite Hn in Eq. ([T]) 

as 

q 

where Iq ~ Z]j e^'^ '^^ Ij and Sq = J2j ^"^'^ '^^^j are the 
Fourier transforms of Ij and Sj , respectively. (From now 
on we set h = 1.) Since A is a small energy scale in our 
case, we can perform a Schrieffer- Wolff (SW) transforma- 
tion in order to eliminate terms linear in A, followed by 
integrating out the electron degrees of freedom. Further- 
more, we can reduce our initial 3D model to a genuine 2D 
problem. The main steps of these calculations are given 
in Appendix[Al We are left with an effective Hamiltonian 
HcS for the nuclear spins in a 2D plane: 

^eff = ^^E^qX"M9) A, (3) 
q 

where 

X„Mq,c.) = -^y^ dte'^-'~^'{[S^{t),S^_^]), (4) 

is a general 2D electron spin susceptibility tensor, 
Xafsiq) = Xa/3(q, w = 0) and g = |q|. TV is the number 
of lattice sites in the 2D plane and a denotes the lattice 
spacing for nuclear spins. Note that (. . .) means aver- 
age over electron degrees of freedom only. We have nor- 
malized X such that it coincides with the density-density 
Lindhard function (see below) in the isotropic and non- 
interacting limit. 

The only assumptions we make are time reversal sym- 
metry of Hq, as well as translational and rotational in- 
variance. The effective Hamiltonian in Eq. ^ is there- 
fore quite general and does not depend on the dimension- 
ality of the system. Note that Eq. ^ is also valid when 
electron-electron interactions are taken into account. It 
is worth emphasizing that the SW transformation ne- 
glects retardation effects. This is appropriate since the 
the nuclear spin dynamics is slow compared to the elec- 
tron one (in terms of energy scales this is related to the 
fact that A <ti Ep). Therefore, electrons see an almost 
static nuclear spin background, and the adiabatic approx- 
imation (for the conduction electrons) is well justified. 
In the case of ferromagnetic semiconductors, such an ap- 
proximation breaks down and retardation effects must be 
taken into account.— If we also assume spin isotropy in 
the 2DEG, then Xap{q,uJ -> 0) = 5a/3Xs(g), where Xsio) 
is the isotropic electron spin susceptibility in the static 
limit. 

In real space, the effective nuclear spin Hamiltonian 
reads 

H.^^-\Y.jfr'lrl'.-. (5) 
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where 



Q/3 

|r| 



-(AV4n,)xa/5(|r|), 



(6) 



is the effective exchange coupUng. The nuclear spins Ir 
are therefore interacting with each other, this interaction 
being mediated by the conduction electrons. This is just 
the standard RKKY interaction^ which, however, as 
we shall see, can be substantially modified by electron- 
electron interactions compared to the free electron case. 



III. ELECTRON SPIN SUSCEPTIBILITY IN A 
2D INTERACTING ELECTRON GAS 

The main result of the previous section is that the mag- 
netic exchange interaction between the nuclear spins is 
mediated by the electron gas. Therefore, the key quantity 
governing the magnetic properties of the nuclear spins 
is the electron spin susceptibility Xs(q) in two dimen- 
sions. The calculation of this quantity in an interacting 
2DEG has been the subject of intense efforts in the last 
decade in connection with non-analyticities in the Fermi 
liquid theory..2126^^^^^^^^ On a more funda- 
mental level, incorporating e-e interactions in the calcu- 
lations of thermodynamic quantities has been an impor- 
tant area of condensed matter theory over the last fifty 
years. In particular, the study of non-analytic behav- 
ior of thermodynamic quantities and susceptibilities in 
electron liquids has attracted recent interest, especially 
in 2D.i2^>2£i21i2ai2ai^i^iHi^iM Of particular importance 
for the following is the recent findings by Chubukov and 
MasloA«2£, namely that the static non-uniform spin sus- 
ceptibility Xsil) depends linearly on the wave vector 
modulus q = |q| for q <s:_ kp in 2D (while it is q^ in 
3D), with kp the Fermi momentum. This non-analyticity 
arises from the long-range correlations between quasi- 
particles mediated by virtual particle-hole pairs, despite 
the fact that e-e interactions was assumed to be short- 
ranged. 

Let us first recall the case of non-interacting electrons. 
In this case, Xs coincides with the usual density-density 
(or Lindhard) response function XLr^ 



XL{q) 



1 



«k,cr — nk+q,o- 



k.cr 



ir] 



(7) 



where rtk is the electron number operator, ek the disper- 
sion relation, and 77 > is an infinitesimal regularization 
parameter. This Lindhard function can be evaluated ex- 
actly and reads in 2D^ 



XL{q)^-N, \l-Q{q-2kF) 



^kl 



(8) 



where Ne = n^jEp is the electron density of states (per 
spin). Note that TVg — m* /it where m* is the effective 
electron mass in a 2DEG. It follows from Eq. ([8]) that 

5xL{q) = XL{q)-XL{0)^Q for g < 2/cf . (9) 



Let us now include electron-electron interactions. It 
is convenient to introduce a relativistic notation with 
p = (po,p) being the (D+l)-momentum where pq de- 
notes the frequency and p the Z?-dimensional wave vec- 
tor (here D — 2). In a zero-temperature formalism, the 
susceptibility can be written diagrammaticallj^ as 



Xs{q) = -jD 2I <^(^'GaiPi-q/2)G^{pi+q/2)Ap,^^,{q), 

pi ,cr,cr' 

(10) 

with a, a' = ±, and where L = aN^'^ is the system 
length, Gaipi) is the exact single-particle Green's func- 
tion and A(g) is the exact vertex function, which can be 
expressed in terms of the exact scattering amplitude T(q) 
as follows!^ 



' E {P2 + q/2)G„, {p2 - q/2) 

P2 



(11) 



This scattering amplitude plays a crucial role as we see 
next. F is for a general scattering event a function of four 
spin variables fii , cr'j , 0-2 , . Nevertheless, one can use a 
convenient parametrization which ensures rotational spin 
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p(Ti(T^cr2Cr2 



,P2 l"'! 0"20'2 ' 

(12) 

where r is a vector whose components are the Pauli ma- 
trices (T^,T!',r^) and r^^^; • r^,^^ = J2 '^a^a'^a^a'- 

a—x^y.z 

Note that F* is spin-independent and corresponds to the 
charge and spin channels, respectively. Following Ref. 
0, we next write the Bethe-Salpeter (BS) equation for 
r~ (corresponding to the spin channel) as follows: 



,M) = {^\„)plP2{<l) + 



(13) 



where iXm)pp'^^) the exact irreducible electron-hole 
scattering amplitude in the spin channel, and 

Rp{q) = -2zGip + q/2)G{p - q/2) (14) 

is the electron-hole bubble4i 

One can exactly solve, at least formally, the BS equa- 
tion ifTS]) using a matrix notation where the matrix in- 
dices run over p. Within this notation i? is a diagonal 
matrix. We find that 

1 



F 



P1P2 / , 
p" 



(^iTT)piP" 



l-r7-Jq)R{q) 



(15) 



P"P2 



This enables us to derive an exact and closed expression 
for the spin susceptibility, given by 



(16) 



PP' 
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In general, T^^ cannot be calculated exactly and some 
approximations are required. The approximation we use 
in the following consists in replacing the exact irreducible 
electron-hole scattering amplitude (r^rr)p,P' aver- 
aged value calculated with respect to all possible values 
of p and p' near the Fermi surface. This is equivalent to 
the following approximation: 



(17) 



We now assume (?o = and suppress the (?o-argument 
in what follows since we are interested in the static prop- 
erties of the spin susceptibility. 



A. Short-ranged interaction 

In this section, we consider a g-independent short- 
ranged interaction potential, which corresponds within 
our notations to rrj,(g) — ~U. This approximation con- 
siderably simplifies the BS equation ifTS)) and the formal 
expression of Xs in Eq. iflGj) . The derivative of Xs{q) 
with respect to q can be expressed in a simple compact 
form: 



1 



dq {l + UIl{q)y 



(18) 



where ^{q) — J2p^pil)/^^- In the low q kp Hmit, 
one can approximate the term Ii{q) in the denominator 
of Eq. ifTSj) by its non-interacting value xl(0) = — A^e- 

The multipHcative factor 1/(1 - UNe f in Eq. ^ 
signals the onset of the ferromagnetic Stoner instability 
when UNf, approaches unity. The Stoner instability is 
supposed to occur for very large ~ 20 according to 
Monte Carlo results.— For smaller Vs < 10, we are still 
far from the Stoner instability. Though this multiplica- 
tive term does not play a significant role at small Ts, it 
increases with Ts, showing the tendency. 



1. Perturbative calculation 

The corrections to the polarization bubble Ii{q) are 
dominated by the first bubble correction to the self- 
energy of Ga-{p)^ These corrections have been calcu- 
lated in second order in U in the small q limit in Ref. [H, 
with the result 



<5n(g) = niq) - n(o) « -Aqxs 



|r(7r) 

37rfct 



(19) 



where xs = lxs(0)| and r(7r) -Um*/4:TT is the (2/cf) 
backscattering amplitude. When UNe <C 1, we recover 
from Eq. (fT8| the known result Sxs{q) — SU{q)^ This 
perturbative calculation therefore gives that Sxs (<?) = ctq 
with a < 0, or equivalently that 1X5(9)! increases with 
q at low g.— Finite temperature calculations along the 
same line imply that Sxs{T) = 7T with 7 < O i^°i^^ This 



result has been confirmed by using a supersymmetric ef- 
fective theory of interacting spin excitations.— 

On the other hand, recent experiments on a low density 
interacting electrons gas in Silicon MOSFETs have found 
that |Xs(2^)l decreases at low temperature^ in apparent 
contradiction with perturbative calculations. This rather 
puzzling situation demands a non-perturbative approach. 



2. Beyond lowest order perturbation theory: 
renormalization effects 

The previous calculations give a spin susceptibility 
which is quadratic in the backscattering amplitude. How- 
ever, it is known that at low enough energy the backscat- 
tering amplitude becomes renormahzed. One should 
therefore instead consider some type of renormalized per- 
turbation theory approximation (RPTA). Shekhter and 
Finkel'stein^ argued recently that the strong renormal- 
ization of the scattering amplitude in the Cooper channel 
may explain the sign of 5xs (T) in the experiment by Prus 
et alM 

Let us introduce T{6,T), the two-particle scattering 
amplitude at a particular scattering angle 9 and tem- 
perature T. The Cooper channel corresponds to 6* = tt 
and we denote by Tc(T) = r(7r,T) the corresponding 
two-particle scattering amplitude. In the Cooper chan- 
nel the two particles that scatter have exactly opposite 
momenta. We first expand Tc in Fourier harmonics such 
that 



rc(r) = E(-i)"rc.n(T), 



(20) 



where n is an integer representing the angular momentum 
quantum. 

In order to analyze the temperature dependence of 
Tc,n{T), it is enough to write the BS equation for Tc 
in a way fully analogous to Eq. lfT3|). In the Cooper 
channel, due to the fact that the momentum of the cen- 
ter of motion of the two scattering particles is zero, the 
integration over the electron-hole bubble gives rise to log- 
arithmic infrared (IR) divergences. These divergences are 
due to particle-hole excitations around the Fermi sea and 
already occur in the ladder approximation {i.e. by con- 
sidering only the bare short-ranged interaction potential 
in the BS equation). We refer the reader to references 
I49I and [13 for more details. The logarithmic IR diver- 
gences can be absorbed by a rescaling of the scattering 
amplitudes in the Cooper channel such that: 



l+Te.,nHEF/kBTy 



(21) 



where Tc^n — Tc,n{^p) is the bare value. The logarithmic 
factors give rise to the well known Cooper instability i^°i^^ 
They are just the 2D equivalent of the one found in the 
discussion of Kohn-Luttinger superconductivity.-^ In Eq. 
(|2T1) , the logarithmic divergence is cut-off by the temper- 
ature at low energy. 
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Another way of recovering Eq. (|2T|) is to regard the 
backscattering ampUtudes F^n as some energy scale de- 
pendent coupHng constants rc,„(A) of the 2DEG (A is 
some running energy scale, the equivalent of T in Eq. 
(|2T1) ). The next step is to write renormalization group 
equations (RG) for these coupHngs. These couplings are 
marginal and the RG equations at lowest order read: 



d\n{K/Ef 



(22) 



The solutions of the RG equations for rc,„(A) are given 
by Eq. (|2T1) where essentially T is replaced by Af^ 

The low energy behavior of the Tc.n{T) strongly de- 
pends on the bare scattering amplitudes being repulsive 
{Tc^n > 0) or attractive (rc,„ < 0). >From Eq. ([211) . 
we immediately infer that when > 0, Tc^niT) 
at low energy/temperature. On the other hand, when 
Ten < 0, Tc,n{T) renormaHzes to the strong coupHng 
regime. Assuming there exists at least one harmonic no 
such that Tc^no < implies that there exists a tempera- 
ture below which the scattering in the Cooper channel is 
entirely dominated by Fc^no-— 

This reasoning relies on the fact that at least one bare 
scattering amplitude in the n*'' harmonic is negative. 
This assumption can be further substantiated with some 
expHcit perturbative calculations of the irreducible scat- 
tering amplitudes j^i^i^ By computing the lowest order 
corrections contributing to the irreducible scattering am- 
plitude Fcirr, (therefore going beyond the ladder approx- 
imation in the BS equation which leads to Eq. one 
can actually prove that there exist higher harmonics such 
that Tc,n < 04^ 

Note that these calculations imply that a supercon- 
ducting instability should develop at a temperature 
T < Tl = £;Fexp(-l/|r,,„J)/fcB . This is en- 
tirely analogous to the Kohn-Luttinger mechanism for 
superconductivity.—"^ Nevertheless, in a typical 2DEG, 
disorder or some other mechanism may provide a natural 
infrared cut-off preventing this superconducting instabil- 
ity to be reached. Let us call A > Tl this infrared cut-off. 
It is worth noting that A should not be too large either 
in order to let Tc^no flow to relative large values of order 
one for this mechanism to be relevant. 

As shown in Ref. [13 this mechanism, if it takes place, 
leads to a non-monotonic behavior of the temperature 
dependence of the electron spin susceptibility and more 
speciflcally to the existence of a scale Tq, below which 
> (note that we use a different sign convention for 
Xs(0) compared to Ref. [13). By fitting the experimen- 
tal data of Prus et al. with such a theory, the estimate 
Tcrio ^ 0.25 — 0.3 was obtained in Ref. [33, which implies 
To ^ 10 K, a surprisingly large temperature scale of the 
order Ep ~ 40 if in this experiment. It is worth em- 
phasizing that such a scale, being dependent on the bare 
value of Fcno, is non-universal and therefore sample de- 
pendent. 

This RPTA also raises a similar issue concerning the 
behavior of the electrostatic spin susceptibility Xs{q) and 



more speciflcally of Tc{q). Following Ref. [43, one can 
write a BS equation at zero temperature for Fc(^), where 
S, is an energy scale in the vicinity of the Fermi energy 
Ep. By linearizing ^ around Ep such that — vpq, one 
immediately infers from l(2T|) that 



, ln{Ep/vpq) 



+ rc.nHkp/2q)' 
(23) 

Note that at flnite temperature and for a realistic system, 
the infrared cut-off A is replaced by max{vpq,kBT, A}. 
This result is also obtained from a renormalization group 
approach by working directly in momentum space at low 
energy where the energy can be linearized around kp. 
At low enough momentum, Tc{0,q) is dominated by the 
harmonic ng such that 



IF, 



1- IF, 



ln(fcF/2g) ■ 



(24) 



By replacing the bare value of the scattering amplitude in 
Eq. (flOl) by the renormalized one, we flnd that Sxs{q) ~ 

— 4gA^e-^^-fer^ also acquires a non-trivial g-dependence. 



In particular, using Eq. l(24|) . we obtain that 



dq 37Tkp^^ ^l^^l^' 



(25) 



which is positive when |Fc(g)| > 1/2. 

Let us assume that A, the aforementioned infrared cut- 
off, is close to the Cooper instability, i.e. ksTj^ ^ A. In 
such a case, the flow of Fc(g) in Eq. l(24|) is cut off by A. 
Using 1/Tc.no = Inl-E-F/fesTi), one obtains that 



Tr ~ : .\ . for q < A/vp 



ln( 



(26) 



is renormalized toward large values. This implies that 
there exists a go ~ lOgL (where qL = kBTL/vp) such 
that j^Sxsiq) > 0, V g < go, and -^Sxs{q) < when 
g ^ go- 

We must point out that the RG equation in Eq. l(25|) 
assumes that all other scattering amplitudes F„ for n ^ 
no can be neglected below some low energy. In order to 
go beyond this approximation, one has to solve (numeri- 
cally) a set of RG differential equations instead. 

This scenario relies on a flne tuning of A compared to 
ksTL (typically this demands Tl ^ A/ks < lOT^). If 
this condition is not met, one may then expect Sxs{q) < 
in accordance with lowest order perturbative calcula- 
tions. Nevertheless, one should mention that an alter- 
native theory, also giving SxsiT) > at low T, has been 
put forward.— This alternative scenario applies for van- 
ishing Cooper amplitudes. In such a case, the anomalous 
temperature dependence of the spin susceptibility is dom- 
inated by non-analytic contributions from particle-hole 
rescattering with small momentum transfer.— Whether 
such a scenario, not considered here, also implies that 
^Xsiq) > is an interesting but open question. 
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The aforementioned considerations immediately raise 
the issue about the ampUtude of qo in a typical inter- 
acting 2DEG. In order to describe the temperature de- 
pendence of Xs(2^), Shekhter and Finkel'stein determined 
Trig such that the experimental behavior of Ref. [H is re- 
produced. This fixes the value of this parameter and also 
the scale ksTo to the order of Ep- One may therefore 
expect qo to be of the order oi kp for a similar 2DEG. 
Another independent way of substantiating these esti- 
mates is to go beyond the ladder approximation in the 
BS equation.— First, as we mentioned before, this al- 
lows one to show that there exists a value no at which 
^c,no < 0. Second, this can give us an estimate for Tl 
and therefore for go- Indeed, in Ref. IH the Cooper 
instability has been estimated to set in at the temper- 
ature ksTL ~ i;i.e-i6^^/('^'«\ where ^ ~ A plays the 
role of the infrared cut-off in Ref. |4^, and the condi- 
tion £,/Ep <^ 1 was assumed. Consistency requires that 
Tl A ^ Ep at large r^, implying qo ~ kp. This is in 
agreement with our previous estimate. 

The previous calculations also give us information 
about the possible shapes for Xs{q)- When renormal- 
ization in the Cooper channel is important, we ob- 
tain at least one extremum around some wave vector 
qo ^ 0{kp). Furthermore, at large q, we should re- 
cover the non-interacting behavior and therefore Xs (<z) ~ 
XL{q) for q oo. Because Xs('Zo) < 0, we expect 
another extremum around a value qi > qo- Since the 
non-interacting behavior is recovered for q >• 2kp, one 
may suspect qi ^ 0(2/cf). From the previous consider- 
ations we therefore conclude that there exist (at least) 
two extrema for the electron spin susceptibility Xs (q) ■ It 
is worth emphasizing that this double-extremum struc- 
ture is a direct consequence of the nontrivial renormal- 
ization of the scattering amplitude in the Cooper chan- 
nel. We have schematically drawn in Fig. [T]the possible 
qualitative shapes denoted by (a) and (b) of Xs{q) — 
Xs{q)/\Xs{0)\ as a function of q/kp and compared it to 
the (normalized) non- interacting XLiq)/Nc at T — 0. 
In the case denoted by (a), we choose Xs{q2) > Xs{0), 
whereas in the case denoted by (b) the absolute value 
of the susceptibility at q2 ^ 2kp is chosen to exceed the 
static value, i.e. Xs{<l2) < Xs(0). The previous consider- 
ations do not allow us to discriminate between these two 
possible shapes of Xsil)- Furthermore, by increasing Ts, 
Xs can evolve from one shape to another. 

On the other hand, if the renormalization in the 
Cooper channel does not take place, e.g. when it is cut 
off by disorder, then the perturbative calculations at low- 
est order apply and give instead that Sxs{q) < at low 
g, 29,30,31 possible shape for Xs{q), consistent with these 
calculations, has been drawn in Fig. [T]and corresponds 
to label (c). We should note, however, that the effect 
of rescattering of a pair of quasiparticles in all different 
channels^ should be carefully examined and may still 
lead to shapes (a) or (b) in Fig. [TJ 



<?() 



T9l 



4 q/kf 




Xs 



FIG. 1: Three possible qualitative shapes denoted by 
(a), (b) and (c) for the normalized susceptibility Xsi<l) = 
Xaiq) /\Xs{0)\ as a function of q/kp (dashed lines) compared 
to the non-interacting value (thin full line). Here, go and 
qi are the positions of the extrema for curves (a) and (b). In 
contrast, the local field factor approximation discussed in Sec. 
IIIIBI results in a monotonic increase of Xs (l) (not sketched 
in the figure), being always larger than the non-interacting 
value. 



B. Long- ranged Coulomb interactions 

In the preceding section, we replaced rj^j.(g) by an al- 
most g— independent constant operator, assuming that 
the Coulomb interaction was screened and, therefore, 
short-ranged. Let us consider in this section the bare 
2D Coulomb interaction, V{q) = jq, where e is the 
electron charge. 



1. Local Field Factor Approximation 

One of the most successful approximations for the cal- 
culation of electron response functions is the local field 
factor approximation (LFFA). It improves the random 
phase approximation^ for which the effective field seen 
by an electron is the field that would be seen by a clas- 
sical test charge embedded in the electron gas. The idea 
of the LFFA to correct the random phase approximation 
and to better account for the correlations existing in the 
electron gas, is to replace the average electrostatic poten- 
tial by a local field effective potential seen by an electron 
with spin a which is part of the 2DEG. (We refer to Ref. 
for a review) . The local field factor G_ can be defined 
as follows: 



G^{q)V{q) = XsHq)-Xl\q), 



(27) 



or, equivalently, the static spin susceptibility Xs can be 
written as 



Xsiq) = 



Xdq) 



l + V{q)G-{q)XL{q) 



(28) 



The precise determination of G_ [q) for all q is an open 
problem. However, the asymptotic regimes, particularly 
the g ^ limit, are quite well estabHshed because they 
are strongly constrained by sum rules4£ In this work, we 
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use a semi-phenomenological interpolation formula given 
in Ref. \M 



G~ (q) « .go 



(29) 



Here, {gfJ-s) "^Xp is the Pauli susceptibility (/^s the Bohr 
magneton and xp > 0); XS = IXs(0)| the renormalized 
value of the spin susceptibility at g = 0, K2 = kprgV^ 
is the 2D Thomas-Fermi wave vector, and go is the pair- 
correlation function at r = 0, describing the probability 
of finding two electrons (of opposite spins) at the same 
position in the electron gas. This phenomenological form 
for G- has been modified from the one originally pro- 
posed by HubbardAS, in order to satisfy exactly the com- 
pressibility sum rule.— >^ The main weakness of this ap- 
proach is the arbitrariness of the chosen form for G_. 
For non-interacting electrons Xp/xs = 1- An approxi- 
mate form for go giving a good agreement with quantum 
Monte Carlo (QMC) calculations has been proposed re- 
cently by Gori-Giorgi et al. and readsi^ 



goirs) « (1 + Ar, + Bri + Cri)e-''^^ /2. 



(30) 



The parameters A = 0.088, B = 0.258, C = 
0.00037, D — 1.46 are fitting parameters reproducing 
QMC results for 50 in a 2DEG.— This approximation 
yields 



Xsiq) 



q + goK2a 

9 + .90^2 (a - 1) 



for q<2qF, (31) 



where a = xs/ixs — Xp) can be regarded as a Fermi 
liquid parameter. The low-g semi-phenomenological ap- 
proximation for the electron spin susceptibility given in 
Eq. jni) results in > 0, Vg, in contrast to the 

lowest order perturbative calculations . 

Note that a direct estimate of G_ (g) by recent QMC in 
a 2DEG gives an almost linear in q behavior up to rather 
large values of g < 2fcF, followed by a more complex non- 
monotonic behavior around 2kF, and finally diverges in 
the large-g limit j^i^i^ This large-g limit is not repro- 
duced by Eq. l(29|) . This is not a serious drawback since 
most quantities of interest are dominated by the low-9 
regime. 

However, the scale q* = goK2ic( — 1) decreases expo- 
nentially with rs according to Eq. l(30|) . This would im- 
ply an almost constant behavior for {q) except at low 
q. When we compare this behavior with available QMC 
data )^^i^^ we find that there is a manifest contradiction. 
Therefore, this raises some doubt about the presence of 
go (a short distance quantity) in Eq. l(29|. 



2. Modified Local Field Factor Approximation 



If we instead replace go in Eq. I|29l) by a parameter 
gi, such that 31^2(0! — 1) = girsy/2{a — l)kF ^ ^kp, we 
have checked that the QMC data for G- (q) are much bet- 
ter reproduced for q's up to 2kp than by the expression 



given in Eq. l(29|) . In such a modified local field factor 
approximation (MLFFA), gi is approximately given by 



51 



Ci 



rs{a - 1) 



Ci 



Xs 
XP 



(32) 



with ^1 some numerical constant of order one. 

We should mention that some other more complicated 
anal ytic al fits of the QMC data have been obtained in 
Ref. [6l|. Nevertheless, we note that the fits used in that 
paper lead to Sxs{q) ~ q^ for 2D, which is in contradic- 
tion with all previous approximations. It seems desirable 
to test Eq. l(32l) with more detailed QMC calculations. 



Comparison of the various approximation 
schemes 



If we summarize the various approximation schemes 
presented in the previous sections, which are perturbative 
or semi-phenomenological, we can clearly ascertain that 
Sxsiq) fx q for q <C qp- Nevertheless, the sign of the 
proportionality constant depends on the approximation 
scheme we used. 

Lowest order perturbation theory in the interaction 
strength leads to dxs{q)/dq < at low g, 29-30 jjowever, 
within the RPTA, renormalization effects in the Cooper 
channel^^^ are important and change the picture given 
by lowest order perturbation theory. In this latter case, 
the RPTA yields an opposite sign for dxs{q)/dq below 
some wave vector qo. The LFFA we used implies that 
dxs{q)/dq > for all q and therefore a monotonic be- 
havior (not shown in Fig. [1]), whereas the RPTA leads 
to a non- monotonic behavior (see Fig. [J). EstabHsh- 
ing a microscopic connection between these two different 
approaches is obviously a rather difficult and open issue. 

The LFFA is a semi-phenomenological approximation 
in which an analytical expression for the local field factor 
G_ {q) is "guessed" with the constraints that the asymp- 
totic behavior should reproduce some known results in- 
ferred from exact sum rules. The unknown parameters of 
G- (q) are fixed from a fit to QMC data.— One may won- 
der whether one can extract some information about the 
possible shapes of Xs{q) directly from the original QMC 
data. The QMC data shows a rather complicated struc- 
ture with two extrema for G-{q) around 2kF (see Ref. 
m or Ref. m (p. 244) ). Though it might be tempting 
to relate the double-extremum structure obtained by the 
RPTA to the QMC results, it turns out to be not possi- 
ble to extract the behavior of Xs{q) from available QMC 
data for G-{q). New QMC calculations directly comput- 
ing Xs{q) instead of G-{q) are thus highly desirable. 

Finally, we have seen that the low q dependence of 
5xs{q) mimics the temperature dependence of 5xs(T) 
which is in agreement with the experiment by Prus et 
al. at low TM- These experimental features may provide 
another, though indirect, consistency check of the RPTA. 
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IV. MAGNETIC PROPERTIES OF THE 
NUCLEAR SPINS 



1. Magnetization and Curie temperature of the nuclear 
spins 



We assume in this section that some nuclear spin order- 
ing actually takes place at low enough temperature and 
analyze how this ordering is destroyed when the temper- 
ature is raised. 



A. Mean field approximation 

Since the interaction between nuclear spins is of RKKY 
type, the interaction is ferromagnetic at short distance 
d <^ kp^ (the large q behavior of Xs{q) is only weakly 
modified by e-e interactions). Furthermore, many mean 
field calculations performed for the 3D Kondo lattice at 
low electron density (neglecting e-e interactions though) 
predict a ferromagnetic ordering . Assuming a low 
temperature ferromagnetic ordering of the nuclear spins 
seems therefore a reasonable assumption. 

We first recall the mean field results for 
completeness.— The Weiss mean field theory gives 
a Curie temperature 



/(/ + 1) 

3k B 4n, 



■XL(g = 0), 



(33) 



where / is the nuclear spin value. 

In 2D this mean field theory yields for T^^^ a depen- 
dence on the ratio ne/rig. For a metal with about one 
conduction electron per nuclear spin, the ratio Ue/ris ~ 1, 
and we recover the result derived more than sixty years 
ago by Frohlich and Nabarro for a 3D bulk metalM 
For a 2D metal the Weiss mean field theory then gives 
ksTc = I{I + 1)A'^/12Ef. For a 2D semiconductor, how- 
ever, the smaller Fermi energy is compensated by the 
smaller ratio rie/ris <C 1. With typical values for GaAs 
heterostructures, / = 3/2, A ~ 90 fieV and a ~ 2 A,— 
we estimate Tc ~ 1 fJ-K, which is very low. (For such 
low Tc's, ignoring nuclear dipole-dipole interactions from 
the start would not be legitimate.) However, this esti- 
mate is just based on the simplest mean field theory and, 
moreover, does not include the effect of e-e interactions. 
It still leads to a finite Tc under which the nuclear spins 
order ferromagnetically. 



B. Spin wave analysis around ferromagnetic 
ground state 



We shall now go beyond the above mean field approxi- 
mation and perform a spin wave analysis. The collective 
low-energy excitations in a ferromagnet are then given 
by spin waves (magnon excitations). 



>From standard spin waves analysis^, the dispersion 
relation of the spin waves in the ferromagnet simply reads 

LUg = I{Jo - J,) = I^a^ixsiq) - XsiO)), (34) 
where Jq is the Fourier transform of Jr defined in Eq. 

At this stage, we already see that the stability of 
the ferromagnetic ground state demands that Sxs{q) = 
Xsiq) — Xs(0) > 0. We can therefore conclude that the 
second order calculation impHes that the ferromagnetic 
ground state is always unstable.™ On the other hand, 
the renormalized perturbation theory approximation de- 
veloped in Sec. IIII A 21 shows that it is necessary to go 
beyond lowest order perturbation theory. 

When renormalization effects are not important in the 
RPTA, the lowest order perturbative results are recov- 
ered, and the ferromagnetic ground state seems unstable 
(though renormalization effects in all channels must be 
carefully taken into account as described in Ref. [33 ). 
When renormalization effects in the Cooper channel are 
important, we expect the two possible shapes denoted 
by (a) and (b) for the static spin susceptibility Xs{q) at 
T = (see Fig. [1]). If the case (b) is favored, then there 
exists a value of q at which ujq < signaling an instability 
of the ferromagnetic ground state. >From this perspec- 
tive, cases (b) and (c) are similar. Another ground state 
must then be assumed and a subsequent analysis is re- 
quired. This will be detailed in Sec. IIV Dl 

On the other hand, if the shape of the susceptibility 
denoted by (a) is favored, the ferromagnetic assumption 
is self-consistent. The RPTA predicts that there exists a 
temperature Tq above which dxs/dq < at low q. This 
implies that there exists another temperature Ti < Tq 
at which the minimum in qi touches the horizontal axis 
signaling an instability. If the Curie temperature Tc is 
larger than Ti, then there exists a temperature regime 
(typically for T > Ti) where the ferromagnetic ground 
state becomes unstable and a different ordering may be 
favored. This case will be analyzed in Sec. IIV Dl On 
the other hand, if the Curie temperature Tc is smaller 
than Ti , the ferromagnetic ground state is self-consistent 
below T^. Such a scenario is in accordance with the one 
obtained from LFFA. Let us therefore analyze this latter 
case. 

The magnetization m per site for a ferromagnet at fi- 
nite T is defined by 

^iT)^I~^^Y.^,^I-^Y.Tl^^^ (35) 



where riq is the magnon occupation number and the sum- 
mation is over the first Brillouin zone of nuclear spins. In 
the continuum Hmit this becomess 

dq 1 



m{T) 



(27r)2 e'^'^^ 



1 



(36) 
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We define the Curie temperature as the temperature 
at which the magnetic order is destroyed by those spin 
waves. This procedure is equivalent to the TyabHkov de- 
coupHng scheme;^ Another way of determining the Curie 
temperature is to analyze at which temperature Tc the 
spin wave analysis breaks down. The Curie temperature 
Tc may be then defined by m{Tc) — 0, which can be 
written as 



1 



a 

T 



dq 



(27r)2 e'^o/'^-s'r<= _ i 



(37) 



For non- interacting electrons in 2D, Xsil) — Xs(0) = 
for q < 2kFr^ where kp is the Fermi wave vector. 
The spin wave analysis, therefore, gives Tc = 0. This 
is in agreement with a recent conjecture extending the 
Mermin- Wagner theorem^^ for RKKY interactions to a 
non-interacting 2D system.— For interacting electrons, 
however, the long range decay of the RKKY interactions 
can be altered substantially and no conclusion can be 
drawn from the Mermin Wagner theorem or its exten- 
sions. 

Let us now include electron-electron interactions (ob- 
tained either by the RPTA or LFFA) . All approximations 
imply that the magnon dispersion is linear in q at low q, 
i.e. that LUq « eg, where 



c = / — 
4 



dxs{.q) 



dq 



(38) 



can be regarded as the spin wave velocity.— Such lin- 
ear spin wave behavior is usually associated with anti- 
ferromagnets while one would expect a quadratic disper- 
sion for ferromagnetically ordered states like considered 
here. This somewhat unexpected Hnear dispersion comes 
purely from electron-electron interactions. 

The perturbative calculations or their extensions to in- 
clude the Cooper pair instability allows us to extract only 
the low q asymptotic behavior of 5xs (?) • Monte Carlo re- 
sults, however, seem to indicate that the local field factor 
G-^{q) is almost linear in q up to g ^ 0[2kF).'^^ We will 
therefore assume that LOq ~ cq for q < q* ^ 0{kF). 

This impHes that for T < T* , where 



T* = cq*/ki 



(39) 



the integral determining m in Eq. I|36|) is entirely domi- 
nated by the Hnear dispersion behavior. Since fast modes 
corresponding to q ^ q* are exponentially suppressed, 
we can easily compute it assuming ujq linear in q for the 
whole q range (extending the upper integration Hmit to 
infinity). We obtain 



m(T) = / [l - {T/Tcf] for T < T*, 



where 



2c 
ksa 



2fc^ 



^ dXsjq) 
im , dq 



(40) 



(41) 



is the Curie temperature. Note that with these estimates 
one has 



T* 
TV. 



aq 



< 1 



(42) 



Such a definition of Tc has been obtained assuming 
that LOq w cq for all q. This approximation has two ma- 
jor aspects: First it regularizes naturally the integral in 
Eq. l(37|) in the UV limit, second only the low energy de- 
pendence ofujq is taken into account, which is consistent 
with a spin wave approximation. 



2. Alternative UV regularization schemes 

In the previous section, we have assumed that ujq « cq 
for all q. On the other hand, one can assume we know 
explicitly LOq for all q in the first Brillouin zone i.e. for 
q < IT ja despite only the asymptotic limits o{5xs{<l) (and 
therefore of ujq) are well established. At large g, we ex- 
pect electron-electron interactions to play a minor role 
and the electron spin susceptibility to be well approxi- 
mated by its non-interacting value, which decreases as 
l/q^ (see Eq. This implies that 5xs{q) ~ -Xs(0) = 

Xs at large q and that the integral in Eqs. l(36|) or (|37|) are 
actually diverging when a ^ 0. If we adopt such a pro- 
cedure, the integral in Eq. (|37|) is fully dominated by the 
short-distance modes, i.e. by the UV cut-off (and there- 
fore independent of any e-e interactions). Such a regu- 
larization scheme is not very satisfying and furthermore 
even inconsistent for a spin-wave approximation which 
relies on the long-ranged modes. Note that the Tc we 
obtain with such procedure is similar (up to a prefactor 
of order unity) to the Curie temperature obtained within 
the mean field theory in Eq. ((SS)) . 

Another regularization scheme consists in cutting off 
the integral in Eq. l(37|) to q < 2C,FkF with Qf a con- 
stant larger than 1. This can be justified by integrat- 
ing out fast modes directly at the Hamiltonian level 
in Eq. ^ since x{l) decreases as l/q^ for q » kp. 
Such a reasoning is equivalent in real space to a deci- 
mation procedure in which a square plaquette containing 
{C,FO'kF)~^ X {(FO-kF)~^ nuclear spins is replaced by an- 
other plaquette with a single average spin. Since at short 
distance, the RKKY interaction is mainly ferromagnetic 
this is equivalent to a mean field procedure. The long 
distance interaction is not substantially modified. The 
main effect of this integration over fast modes is that the 
UV cut-off in Eq. (|37|) is now of order 2kF instead of 
ir/a. Although such a procedure does not allow for an 
exact calculation of Tc since Uq is not known around 2kF, 
this considerably boost the Curie temperature by orders 
of magnitudes compared to the previous regularization 
scheme and in the same range as the Curie temperature 
determined from Eq. (|4T|) . 

In the following, we will therefore use Eq. l(4T|) as our 
definition of Tc. This has the advantage of providing us 
with a simple closed formula. Furthermore, such a Tc 
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is consistent with the long range approximation for spin 
waves. 



Self-consistent calculation of the Curie 
temperature 



3. Numerical estimate of the Curie temperature 



Eq. (|4T|) gives us an estimate of the Curie temperature 
as a function of the derivative of the electron spin sus- 
ceptibility. We have computed these quantities in Sec. 
mil for various approximate schemes. We are therefore 
now ready to give estimates for the Curie temperature 
Tc. Let us start with the RPTA. Assuming Tc ~ 0(1) 
in Eq. l(26|) . we obtain a Curie temperature ~ 20 fiK 
for typical GaAs 2DEG parameters with ~ 1. Setting 
q* ~ 2kF in Eq. I|42p and using the same parameters, 
one obtains T* Tc/50. 

At larger r^, Tc is enhanced for two reasons: first kp^ 
increases linearly with r^, and second, the value of the 
spin susceptibility at q = 0, xs — IXs(0)|, which is es- 
sentially the Pauli susceptibility at small r^, increases 
linearly with r^.--- An approximate value of xs can be 
extracted from QMC calculations.— One obtains, for ex- 
ample, Tc 0.3 mK for = 5 and Tc ~ 0.7 mK for 
Ts = S. One may even obtain larger values of Tc in the 
mK range for larger Vs since Tc increases quadratically 
with Tg. Furthermore, when UNc is no longer negligible 
compared to 1, Tc (resp. T*) is even further enhanced by 
an additional factor 1/(1- C/7Ve)2 (see Eq. (UHl). Close to 
the ferromagnetic Stoner instability of the electron sys- 
tem, reached when UNc '■-^ 1, the Curie temperature Tc 
for the nuclear system is strongly enhanced as could have 
been anticipated. 

On the other hand, if we use the local field factor ap- 
proximations developed in Sec. IIII B[ we can determine 
Tc by inserting Xsil) obtained from the LFFA into Eq. 
(SIl): 



T = 



I A 31 



2kf 



(43) 



where gi = go [Eq. ^] or g^ = gi [Eq. The en- 

ergy scale (a — l)'^giV{a) can be interpreted as a renor- 
malized screened potential due to collective interaction 
effects that are incorporated in the LFFA. 

If we use gi = go, the LFFA gives an exponential en- 
hancement of Tc with increasing the interaction parame- 
ter rg. Yet, as we have discussed in Sec. IIII B| this form 
for the local field factor G_ cannot really be trusted at 
large > 1 when compared to QMC data. If we instead 
use the MLFFA, which seems to be in better agreement 
with QMC data, then gi = gi and we obtain Tc ~ 0.4 mK 
for Ts ~ 5 and Tc ^ 1 mK for ~ 8. These values are 
consistent with the ones found using renormalized per- 
turbation theory. 

Note that the ratio A/ (a — l)^goV{a) can be regarded 
as the small parameter of our theory. For some large 
value of Ts, the dimensionless parameter A/ {a — l)'^gV{a) 
may approach unity and the truncation of the Schrieffer- 
Wolff transformation at lowest order becomes unjustified. 



If we assume a ferromagnetic ordering of the nuclear 
spins, this generates some rather large nuclear Over- 
hauser magnetic field Bes{m) ~ 0(1 T) in a GaAs 2DEG 
for m < / that we have neglected before. Therefore, the 
spin degeneracy of the conduction electrons will be lifted 
by the Zeeman energy g/isBeff which should have some 
effects on the determination of Tc, too. A self-consistent 
procedure including feedback effects is therefore required. 
This is the purpose of this section. 

The electron spin susceptibility is no longer isotropic. 
In the following we will assume that it is still diagonal but 
has now a longitudinal component xt (in the direction of 
the magnetic field produced by the nuclear spins) and 
a transverse component x^. Taking into account such 
an anisotropy, our SW transformation is still valid and 
we obtain an effective spin anisotropic 2D long-ranged 
Hamiltonian, which replaces Eq. ([5]): 



JZ TZ TZ 

''|r-r'|-'r ^r' 



-'ir-r'l 



(7+4^+/-/?) 



where, in q space, Jq 



z/± 



(44) 

-(^V4ns)x^/^(g). The dis- 



persion relation for magnons becomes now 



(45) 



and therefore generically acquires a gap for Bcs ^ 0. Let 
us denote by A(_Bcff ) = A(m) such a gap. The expansion 
of the magnon dispersion for small q leads to 



I A 



3 

^9 



I A 



cq, 



(46) 



where we have assumed a linear in q behavior for x't (<z) 

and defined c' — I-g^{q = 0). We assume ferromagnetic 
order and so c' > 0. In the following we develop a simple 
approach of the Landau type. Eq. ((37|) which determines 
the Curie temperature is still valid and now reads: 



1 = 



a 

T 



(27r)2 e/5^^+/3c'|q| _ i' 



(47) 



where [3 — l/ksT. As before, we assume the linear dis- 
persion to be valid for all q, which naturally regularizes 
the integral in the UV limit. The latter integral equation 
can be rewritten in a more compact form as 



1 



27rJ {Pd) 



rLi2(e" 



(48) 



where Li2 is the dilogarithm function. Let us first focus 
on the magnetization m. >Prom Eq. l(36l) we have 



m(T) 



1 



27r (/3c') 



rLi2(e" 



-/3/A(m(T))^ 



(49) 
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which gives an electronic Zeeman energy scale A(m = 1) 
of the order I A ~ 0(1 K). This scale is much larger 
than the typical Curie temperature we estimated before. 
In the ferromagnetic phase, we are therefore in a regime 
in which A(m = 1) ^ fcsT, where T is the temperature. 

In this regime, one can use Ref. [s^ to estimate A' ~ 
A^/Ep{akp)'^rl, where Tc is the renormaHzed scattering 
amplitude defined in (|24)l . Using Eq. l(37|) for the Curie 
temperature Tc, we can estimate b ~ {akp){A/ Ep) ^ 1. 
Our problem thus corresponds to a regime in which feed- 
back effects are negligible, which justifies the previous 
isotropic approximation for the electron spin susceptibil- 
ity. 



FIG. 2: Magnetization as a function of to = (T/T^) - 1 for 
various values of the dimensionless parameter b. For small 
values of b, the behavior described by Eq. l|40|l is recovered. 



When T Tc, to ^ and hence A(to) 0. Let us 
denote by r° the Curie temperature obtained from Eq. 
(|37|) by neglecting Bcs. We assume that the temperature 
is below but close to Tc. We can therefore expand the 
gap as A = A'm + • • • with A' = ^^(m = 0). The 
dilogarithm around A = then expands to: 



Li2(e 



6 4 



72 



f\3 



(50) 



where g{x) = a;(ln2; — 1). Introducing the dimensionless 
parameters t = T/Tc - 1, io = T/T° -1,6 = lA'/ksTc, 
the self-consistent equation (|49| for m{T) becomes 



{t) = -2to + -^bm {1 ~ \n{bm)) + 



3P 
27r2 



127r2 



(51) 



This equation can be easily solved numerically. Note 
that by integrating this equation with respect to to, one 
obtains the Landau functional. 

In Fig. [21 we plotted the magnetization as a function 
of to for various values of the dimensionless parameters 
b. Two interesting features can be noticed: First, Tc 
increases compared to T^ by increasing b. Second, the 
magnetization rises sharply as T is lowered through T^ 
which is reminiscent of a first order phase transition for a 
system in the thermodynamic limit. This effect becomes 
pronounced at values of 6 ~ 0.15. 

Let us now estimate the value of b in our system. To 
this end, we need to estimate A' which can be regarded 
as the susceptibility of the susceptibility. Since the limits 
T ^ and Soft — ^ for the susceptibility Xs{T,Bcs) 
do not commute,— a little care is required. When the 
nuclear spins in the GaAs heterostructure are polarized, 
they generate a rather large effective magnetic field Bcf / 



D. Spin wave analysis around a helical ground state 

We have assumed so far that the ground state is fer- 
romagnetic. This assumption obviously depends on the 
sign of Sxs{q)- We have seen that the RPTA predicts 
two possible scenarios: Either 5xs{q) < in some range 
of g (case (b) in Fig. [T| or 5xs{q) becomes negative above 
some temperature Ti. A positive sign, on the other hand, 
is always obtained from LFFA. In contrast, second order 
perturbation theory gives a negative sign for 15x8(9), in- 
dependent of the magnitude of q. If we assume that there 
is a temperature range where 6xs{q) < 0, this clearly im- 
plies that the assumption of ferromagnetic order is invalid 
since the magnon spectrum has ujq < 0. This corresponds 
to the situations (b) or (c) in Fig. [H which may depend 
on the Fermi Hquid parameters or arise when T exceeds 
some temperature Ti. 



1. A new ground state 

A different ground state for the nuclear spins thus has 
to be assumed. If the minimum of Sxs{q) is reached at 
some finite wave vector qi ^ kp, a reasonable ground 
state of such a system could be an incommensurable he- 
lical state, or a spiral state, with a wavelength of the helix 
of the order of Ai = 27r/gi. 

We stress that, on the scale of the lattice spacing, the 
RKKY interaction remains ferromagnetic, and the neigh- 
boring nuclear spins are still ferromagnetically aligned. 
On the scale of the Fermi wavelength, however, opposite 
spin alignment is favored. A gain in energy may be ob- 
tained by a slow variation of the nuclear spin alignment, 
i.e. by helical order. If such a nuclear spin arrangement 
is reached for some temperature range, it implies that the 
total magnetization vanishes over distances much longer 
than Xp. 

We assume that the ground state of the nuclear spins 
is described by a local magnetization rrii = (1^) that is 
confined to the spin {x,y) plane and that rotates in spin 
space around the z-direction. A full rotation is described 
by the wavevector qi . 
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These rrii play the role of the order parameter in the 
system. In order to compare with a ferromagnet, we can 
characterize the order by locally rotating the mi so that 
they all map onto the same mo: 



nuclear spins can be written as 



1 J 
Has = --J2ill,l!j!)R.3.,B^ I] 



(53) 



m.i = i?i mo = 



cos(ri • qi) sin(rj • qi) 0^ 
-sin(rj • qi) cos(rj • qi) I mo, 











1, 



(52) 

where mo = m(ri ~ 0) is the magnetization at site 
Ti = 0. Note that mo has the same interpretation as 
the average magnetization of a ferromagnet. 

The helical ground state breaks the translational as 
well as the rotational symmetry of the Hamiltonian in Eq. 
([3]). In particular, the ground state is degenerate with 
respect to a change of the direction of qi in the 2DEG 
plane. With the choice of a specific direction qi, this 
symmetry is broken. There exists, therefore, a Goldstone 
mode that tends to restore this rotational symmetry by 
inducing changes on the direction of qi . Obviously, such 
a mode would destroy the assumed order immediately. 
Realistically, however, the specific helical state can be 
pinned by the system, such as by a disorder lattice con- 
figuration or by some Dzyaloshinskii-Moryia interactions, 
and the Goldstone mode becomes massive.— In the fol- 
lowing we will assume that this is the case. Concretely 
this means that excitations with small q in the perpen- 
dicular direction of qi are strongly suppressed. 



In this new basis, the spin-wave analysis is analogous to 
the ferromagnetic case and rather standard (see e.g. Ref. 

The ground state energy of the helimagnet then be- 
comes 



TP JXX 



(54) 



Let us compare Eq. (|54|) to the ground state energy of a 
ferromagnet. If all spins are aligned along the x direction 
the ground state energy of the ferromagnet is 



TpFM 



7^2 

_ \ ^ 7^a: 
2 ^ 



jXX 



(55) 



The energy of the helical state is thus lower than that of 
the ferromagnet if J^^ > J^^q. Furthermore, the helical 
state has the lowest energy at the wave vector qi , where 



J^^ — jvy has its maximum. 



The low energy excitations above the ground state can 
also be obtained in a straightforward manner from Eq. 
([53)) . We find the following two branches of the spin wave 
spectrum. 



2. Spin wave analysis 



/I) = 
(2) _ 



(J, 



JXX \ 

Qi "^Iqi+qK' 



f JXX J^^\ 

Wqi '^Q )■ 



(56) 
(57) 



Our starting point is the effective spin Hamiltonian 



Q/3 



Clearly, Wq^^ = at q = 0. On the other hand, we 
see that the helical ground state can only be stable if 



,(2) 



of Eq. ^ with an exchange interaction matrix 
j"/3(|r^ - r^l) defined in Eq. © that can be off-diagonal 

in general. Since J"^^ is proportional to the spin suscepti- vicinity of q 

bility, as can be seen from Eq. (|A6p . we can assume, how- case. The second branch, ujq \ then has a gap. 



> 0, which means that must not exceed J^^ in 
0. If J^^ = J^^ this is indeed the 



ever, that Jfj' 



J^j and that Jj"'' = — j/^" iov a ^ [3. 
These assumptions are related to the conservation of the 
total spin of the electron system?^ 

A convenient way to perform the spin wave analysis 
is to first perform a local rotation of each nuclear spins 
Ij as in Eq. (|52l) such that they all become parallel to 
each other like in a ferromagnet. We thus define a local 
(right-handed) set of axes described by the unit vectors 
e^,ef ,ef, with being parallel to mi.-- In principle, ef 
and e? can be chosen arbitrarily. It is convenient though 
to choose parallel to the spin rotation axis z. Then we 
can write li = llej+lfef + lfe^. These new components 
/f are connected to the original components /f through 
the matrices Ri as {il,if,if)'^ = R\{If , If , I^Y , where 
T denotes here the transposition. Let Jy be the 3x3 
matrix associated with JfJ^ , then the Hamiltonian for the 



3. Effect of gapless modes 

We see that there is a gapless mode Wq in the system, 
given by Eq. ijSG]). If q is such that |qi -|-q| = qi, then ujq 
remains strictly zero. This is the aforementioned Gold- 
stone mode. Such fiuctuations, therefore, are assumed to 
acquire a mass, associated with an energy scale Ag. In 
the longitudinal direction, however, where q is parallel 
to qi, Wq grows for |qi -l- q| > gi and |qi -t- q| < gi pro- 
portionally to the increase of Xs(|qi with respect to 
its minimum. 

The minima of the spin wave spectrum are pushed to 
finite momenta (?, notably to the minimum of the spin 
susceptibility Xsil) ^it finite q = qi- Provided Xs{q) re- 
mains analytic around the minimum at qi, it no longer 
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grows linearly but, in general, as {q — qi)"^. The system, 
therefore, can no longer benefit from a linear Xs{q) to sta- 
bilize the ferromagnetic order and different, nonuniversal, 
energy and length scales will affect the Curie temperature 
Tc (see below). 

We have seen above that the quantity mo — |mo| can 
be used as the order parameter of the heHcal state. Every 
spin wave reduces toq , and we can thus use again Eq. l(36|) 
in order to express how the average local magnetization 
is reduced by the spin excitations. 



niQ = I — a 



dq 



(27r)2 - 1' 



(58) 



For q ^ 0, the integrand becomes singular as l/q'^M 
The singularity cannot be compensated by the factor q 
of the spherical integration measure q dq as in the pre- 
viously discussed ferromagnetic case. This is the same 
situation that is met when calculating spin wave ex- 
citations for systems without long-ranged interactions. 
There the singularity is directly linked to the Mermin- 
Wagner theorem.— In the present case, however, the 
small q values are cut off at some finite inverse length 
scale tt/L, associated with the energy A, the infrared 
cut-off frequency introduced in Sec. IIIII 

The scale L lifts the singularity in the integral l(58l) by 
cutting off the momentum at q « tt/L, and so by effec- 
tively opening a gap for the excitations. The spin system 
can maintain a quasi-order over the length L. An ad- 
ditional reduction of the singularity also arises from the 
fact that the system is not truly 2D but a layer with a 
finite width Wz, containing about 50-100 planes of nu- 
clear spins (see also App. A). This length scale, however, 
must be compared with the typically much longer wave- 
lengths of the spin waves, and thus can only account for 
a partial regularization. The singularity is dominated 
by the minimal curvature of ujq in all directions of q. 
In the present case there are two main directions, the 
longitudinal one parallel to qi, where Wq — Cq^, with 
C = —d'^Jq/dq^\q=q^, and the transverse one perpendic- 
ular to qi, where the curvature is imposed by Aq. For 
the stability of the ground state we must assume that 
the pinning strength Aq is large compared to the ener- 
gies imposed by the J^^ . The singularity in Eq. (|58l) is 
thus dominated by 1/ (iCq^ 

Let us assume L ~ 10/im. Cutting off the upper inte- 
gration limit hy kp^ the singular part of Eq. l(58|) becomes 



77^ ln(fci?L/7r). 



Cpq CP 



(59) 



The logarithm yields a factor exceeding one. The heli 
cal order cannot be stable if the expression in Eq. l(59] 
becomes larger than one. This allows us to define a tem- 
perature 



C 



la? InikpL/TT) 



(60) 



above which the gapless mode definitely destroys the he- 
lical order. 

We see that much of the stability depends on the 
value of C, which means, on the curvature of Xsiq) 
around its minimum at qi. We can very roughly esti- 
mate C/a? ^ Jq^KakpY , and see that this temperature 
T'-^ can actually be quite high. Let us now write the 
integral in Eq. (|58l) in the form 



dq 



InikpL /-k) 

Jc 



dq 



|q|>/CF 



(27r)2 eP'^'i - 1 ■ 



(61) 



Let us now further introduce a temperature T**, similar 
to the temperature T* for the ferromagnetic case, below 
which the integral determining mo is entirely dominated 
by the quadratic dispersion behavior. This means we 
assume that cjq « Cq^ (in the direction parallel to qi) 
up to a q** < qi, kp. We can set 



T* 



C{q**f/ki 



(62) 



For Ti < T < T**, therefore, Eq. ^ is controlled 
entirely by the first logarithmic part. We thus obtain 



^„«/_„2MfcFi/-).„^_.,. ^^^^^ 



-fcsT = /(I - T/T'^ 



(63) 



valid for T < T** . 

The temperature T'~^ can be seen as a generalization of 
Tc for the helical case. It differs from the ferromagnetic 
case through its dependence on external, non-universal 
cut-off scales. This loss of universality is an essential dif- 
ference compared to the previously studied ferromagnetic 
order. This also indicates that the helical order is much 
more fragile with respect to external conditions than the 
ferromagnetic order. 

At higher temperatures the formation of defects and 
magnetic domains will further tend to destabilize the or- 
der as well. It is therefore possible that the helical order 
is destroyed well below T'^. 



V. CONCLUSION AND DISCUSSION 

Summary In this paper we have examined the interplay 
between an interacting electron liquid in 2D with mag- 
netic order in a lattice of nuclear spins. We have based 
our investigation on a Kondo lattice model, in which the 
electrons couple weakly to the nuclear spins through the 
hyperfine interactions. In this way, an effective coupling 
of the RKKY type is induced between the nuclear spins 
which, as obtained through a Schrieffer- Wolff transfor- 
mation, is expressed in terms of the static electron spin 
susceptibility Xsil)- 

Electron-electron interactions in 2D can substantially 
modify the shape of Xs (?) and therefore profoundly affect 
the magnetic properties of the nuclear spin system. A 
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magnetic order can arise because the conditions for the 
Mermin- Wagner theorem are not met due to the long- 
range character of Xs (r) ■ Much depends thus on the pre- 
cise shape of Xsi<l)- Based on a renormaHzed perturba- 
tion theory we argued that for short-ranged interactions, 
Xs (q) should have the forms sketched in Fig. [TJ When 
renormalization of the scattering amplitudes is impor- 
tant, Xs{l) has two extrema at values qo,qi '■-^ kp, which 
lead to the two generic situations labeled by (a) and (b) 
in Fig. dJ If renormalization effects are unimportant, 
this leads to the form (c) as sketched in Fig. [TJ The dis- 
tinction between (a), (b), and (c) is non-universal and, 
presumably, depends on the sample disorder, the inter- 
actions, and the temperature. 

For long-ranged Coulomb interactions, on the other 
hand, a calculation based on a local field factor ap- 
proximation produces a monotonic increase of Xs (<?) (not 
sketched in Fig. [1]). 

Such a monotonic increase or the case (a) stabilize a 
nuclear ferromagnet. In the case (b) the ferromagnetic 
order has an instability at the wave vector qi ~ kp, which 
corresponds to the absolute minimum of Xs((/), but fer- 
romagnetic coupling is maintained at short distances (as 
compared to the Fermi wavelength). Similar behavior 
emerges also for case (c). We argued that the nuclear 
ground state then has a quasi-order which is (nearly) 
ferromagnetic on short distances but rotates the local 
magnetization on a scale of 1 /qi , thus providing a helical 
order. 

Experimental implications The transition temperature 
Tc describes the temperature above which the magnetic 
order breaks down. While the mean field approximation 
predicts a very low transition temperature of the order of 
~ l^K, we have seen that electron-electron interactions 
can considerably increase the value of Tc- From our var- 
ious approximation schemes we obtain a Tc in the rriK 
range for both short-ranged and long-ranged interactions 
for = 5 - 10 (see Sec. IIVB3L 

If temperature is decreased below Tc, the nuclear spins 
order and generate an effective magnetic field BcS, which 
can be very large in GaAs 2DEGs, this being in con- 
trast to Si MOSFETs, which have a much smaller Over- 
hauser field. This internal magnetic field has impor- 
tant consequences for the thermodynamic behavior of the 
electron spin susceptibility, which can be derived from 
the non-analytic dependence of Xs{T, Beg)on temper- 
ature and magnetic field following Ref. [SJ (the lim- 
its T and BcS ^ do not commute). In a Si 
MOSFET 2DEG, we expect gfiBBcS < ksTc and there- 
fore the calculated and observed^^ linear in T behav- 
ior of Xs[T) is also valid below Tc. However, in GaAs 
2DEGs one has clearly gfisBcS S> ksTc, and one may 
therefore expect (xsiT) — Xs(0)) oc T above Tc and 
{Xs{T) -Xs{0)) cx m cx [1-{T/Tc)^] below Tc in contrast 
to Si MOSFETs. This implies an upturn of Xs(T) around 
T = Tc in a GaAs 2DEG. 

Open questions There remain many open questions. 
Mainly, a detailed study of Xs(q) at values q 2kp 



could clarify the scenarios qualitatively sketched in Fig.[TJ 
Monte Carlo simulations and experiments may provide 
further insight here. It would be desirable to establish a 
general magnetic phase diagram for the nuclear spins as 
a function of and T . Possible new phases such as a 
nuclear spin glass phase are likely due to the complexity 
and richness of the problem. Finally, disorder may play 
an important role by providing a further cut-off length, 
and the interplay with electron-electron interactions re- 
quires a separate investigation. 



Acknowledgments 

We would like to acknowledge fruitful discussions with 
A. Millis, G. Vignale, and useful correspondence with A. 
Finkel'stein. This work is supported by the Swiss NSF, 
NCCR Nanoscience, and JST ICORR 



APPENDIX A: DERIVATION OF THE 
EFFECTIVE MAGNETIC LOW ENERGY 
HAMILTONIAN 

1. Effective Hamiltonian 

We start from the general Kondo lattice Hamiltonian 
in Eq. ^ with Hdd = 0. Since A is a small energy 
scale in our case, we can perform a Schrieffer- Wolff (SW) 
transformation in order to eliminate terms linear in A, 
followed by "integrating out" the electronic degrees of 
freedom. Keeping the lowest order terms in of the 
SW transformation, we are left with an effective interac- 
tion Hamiltonian Hsw that reads 

Hsw =Ha-\[S,[S,HQ\]. (Al) 

S is defined by Hn + [S, Hq] — 0, which is solved as 
S — LQ^Hn, where Lq is the Liouvillian superoperator. 
Let us define 

U^^[S,[S,Ho]], (A2) 
which can be rewritten as U — ^[Lq , iJ„]. Using an 

oo 

integral representation for Lq ^ = / die*'^" one 



obtains 

C/ = -^^ dte-^'[Hn{t),H,,], (A3) 

where H,,{t) = e'^"*H„ = e'"°^ H„e-'"°\ and r? -> 0+ 
ensures the convergence of the time integration. Using 
the definition of Hn given in Eq. ((2]), we are left in q- 



16 



space with 

^ q,q' ^ 

■ yi2 fOO 

= -hpT. / dte-^'{l^4iS^(.t)^<'] (A4) 

' q,q' " 

Summation over Greek indices is implied. The first com- 
mutator enters the definition of the susceptibility in Eq. 
(|A6P below. The second commutator can be straightfor- 
wardly computed by going back to real space. We obtain 
= i^^^^Iq+q' where e"^'^ is the fully antisym- 
metric tensor. 

We next take the equilibrium expectation value over 
electronic degrees of freedom only, denoted by (...). 
Furthermore, we assume translational invariance in the 
2DEG, which implies (OqOq') = A^(OqO_q>(5q+q',o, with 
N being the number of sites in a 2D lattice. Since Hq has 
time-reversal symmetry the term proportional to e"^T/^ 
in Eq. l|A4p drops out. Together with the reduction to 
a 2D problem (discussed below), this allows us to bring 
(U) to a much simpler form, 

iU) = ^a'J2^^xMq)I^-.,, (A5) 
q 

where Us — is the 2D nuclear spin density. The 
quantity 

/"CXD 

X«M9,^) = -^y^ dte-^-'-^mS^{t),S^_^]) (A6) 

is the 2D electron spin susceptibility which includes 
electron-electron interactions. 



2. Reduction to 2D 

As already mentioned, the electron gas is not strictly 
2D but has a finite thickness, typically of the order of 
Wz ^ 5 nm and therefore contains several layers of 2D 
nuclear spin planes. However, using that the electron 
wave function is confined in the third dimension {z- 
direction), we can reduce the original 3D system to an 
effective 2D one. Let us suppose a lattice site is labeled 
by Tj = {rj//,Zj), where (r^// is the planar coordinate 



and Zj the position in the perpendicular direction. We 
consider ~ 50 layers such that Ni = N x N^. 

If we assume that Wz is sufficiently small so that the 
electrons are confined in a single mode (j>{z) in the z- 
direction, all nuclear spins in a column along the z- 
direction couple to the same electron wave function. 
Fluctuations of the nuclear spins along this direction are 
expected to be weak. In a mean-field-like description, we 
may thus replace the column by a single 2D nuclear spin 
as follows. Let us separate out the perpendicular mode 
from the electron spin operator as 

S,(r,//,z,)-^^M^S,(r,//), (A7) 

Wz 

with Si{ri//) being a 2D electron spin operator. The 
mode ipizi) can then be used as an envelope function for 
the nuclear spins 

I»(rV/) = — / dz^\(P{z^)\%{r,//,Zi). (A8) 

Wz J 

Since Wz is determined by Wz — J dz\(j>{z)\'^ , these new 
operators Ii(rj//) satisfy the standard spin commutation 
relations. The remaining Hamiltonian is now strictly 2D. 

Alternatively, we can argue as follows. Since kp^'3> 
Wz, we can approximate jr^ — rj| w jr^// — rj//| and 
therefore 

xj-'^dr^/z-r,//!) |^^^//(r,.//,z,.)j ■ (A9) 

This amounts to replacing a single nuclear spin 
at position (jj^^^Zj) by an average nuclear spin 
-7= -^j)- This is fully consistent with the fact 

that the RKKY interaction is almost constant and ferro- 
magnetic in the z-direction perpendicular to the 2DEG 
(this is at least the case for the RKKY interaction ob- 
tained by neglecting electron-electron interactions but we 
expect that the interactions do not modify significantly 
the RKKY interactions in the z-direction). 

Our problem has now been reduced to a 2D system con- 
sisting of TV = Ni/Nz nuclear spins interacting with long- 
ranged interactions. The effective nuclear spin Hamilto- 
nian Hcs = {U) is finally given by Eq. ([3]). 
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